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Abstract 

Molybdenum (VI) biosorption process was investigated by marine algae Cystoseria indica pretreated with 0.1 M 
CaCb solution in a packed bed column. The biosorbent was characterized by FTIR, BET and SEM analyses. The 
results showed that Mo (VI) ions should be chelated with the hydroxyl, carboxyl and amine groups of the biomass. 
The effects of inlet metal concentration and flow rate on biosorption process were investigated and the 
experimental breakthrough curves were obtained. Results showed that the maximum biosorption capacity of 
Ca-pretreated C. indica for Mo (VI) was found to be 18.32 mg/g at optimum flow rate of (1.4 mL/min). The 
controlled-rate step shifted from external to internal mass transfer limitations, as the flow rate increased. Also, it was 
observed that the breakthrough and exhaustion time decreased from 17.14 hr to 9.05 hr and from 0.006 h to 
0.002 hr respectively, with the increase of flow rate from 0.7 to 2.1 ML/min. The increase in the initial concentration 
of Mo (VI) solution from 30 to 95 ml min"^ increases the adsorption capacity from 18.32 to 30.19 mg/g and 
decreases the percentage of Mo (VI) removal from 61 to 38%. Also, the treated volume was the greatest (1.42 L) at 
the lowest inlet concentration. Column data obtained under different conditions were described using the Thomas, 
Yoon and Nelson, Yan and Belter models. The breakthrough curve predictions by Belter model were found to be 
very satisfactory. 
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Introduction 

Continuous discharge of industrial, domestic and agricul- 
tural wastes in rivers and lakes causes deposit of pollutants 
in sediments. Such pollutants include heavy metals, which 
endanger public health after being incorporated in food 
chain. Heavy metals are present in different types of 
industrial effluents, being responsible for environmental 
pollution [1]. Incidence of heavy metal accumulation in 
fish, oysters, mussels, sediments and other components of 
aquatic ecosystems have been reported from all over the 
world [2,3]. 

Many types of yeast, fungi, algae, bacteria and some 
aquatic plants have been reported to have the capacity 
to concentrate metals from dilute aqueous solutions and 
to accumulate them inside the cell structure [4-6]. 
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The potential application of microorganisms as biosor- 
bent for the removal of heavy metals has been recognized 
as an alternative to the existing conventional methods for 
detoxification of industrial wastewaters [7,8]. 

Molybdenum is an essential trace element for both 
plants and animals [9]. Molybdenum deficiency has often 
been reported but in large concentrations, molybdenum 
may be toxic as it leads to secondary copper deficiency 
[10]. Therefore, the removal of this ion from wastewater 
and groundwater is of significant importance from an 
environmental point of view. 

One of the heavy metal removal technologies called 
biosorption uses non-living microorganisms in treating 
wastewater containing heavy metals [11,12]. 

The role of various microorganisms by biosorption in 
the removal and recovery of heavy metal(s) has been well 
reviewed and documented [13-18]. It is particularly the 
cell wall of certain algae, fungi and bacteria which is found 
to be responsible for the phenomenon of biosorption [19]. 
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Among these, the brown algal biomass has been shown 
to be highly effective, reliable, and predictable in the re- 
moval of heavy metals from the aqueous solutions [20-22]. 

The packed bed adsorption has a number of advan- 
tages: it is simple to operate, attains a high yield and it 
can be easily scaled up from a laboratory-scale proced- 
ure. The stages in the separation protocol can also be 
automated and high degrees of purification can often be 
achieved in a single step process [23]. 

In our previous study [24], we showed that the pretreat- 
ment of Cystoseria indica by 0.1 molar CaCl2 solutions 
improved the uptake capacity of uranium more than 30%. 
In the present work, the molybdenum (VI) biosorption by 
the Ca-pretreated brown algae {Cystoseria indica) was 
investigated in a packed bed column. The effect of design 
parameters, such as flow rate and influent concentration 
on molybdenum (VI) biosorption were examined. For a 
proper design of a biosorption column, an accurate pre- 
diction of the breakthrough curve is needed. Therefore, 
the experimental results obtained from the continuous 
flow packed bed column were fitted to the Thomas, Yoon 
and Nelson, Yan and Belter models. 

Materials and methods 

Preparation of biosorbent 

The brown alga, Cystoseria indica was obtained from 
The Persian Gulf on the coast of Qeshm, It was thor- 
oughly washed with distilled water and sun-dried on the 
beach and in an oven at 50°C overnight. The dried bio- 
mass was ground in a laboratory blender and sorted 
using the standard test sieves. The batch of biomass with 
particle size of 1.0-1.25 mm was selected for subsequent 
Ca-pretreatment. Ca-pretreatment of the biomass was 
carried out as follows: a sample of 10 g of biomass was 
treated with 0.1 M CaCl2 solution (1000 mL) for 
12 hours under slow stirring. Then the biomass was 
washed with distilled water. Finally, the biosorbent was 
again dried at 70°C for 8 h. 

Chemicals 

The standard solutions of Mo (VI) were prepared by di- 
luting 1000 mg/L Mo (VI) solution by dissolving analyt- 
ical grade sodium molybdate (Na2Mo04.2H20) in water. 
Diluted solutions were prepared at room temperature in 
distilled and deionized water. The pH of influent solu- 
tions was adjusted with a pH meter (Metrohm, Model 
780) by using 0.1 M HCl and/or 0.1 M NaOH. 

Continuous flow column system 

All of the experiments were performed in a packed-bed 
column at room temperature (25 ± 2°C) and pH = 2. A 
preliminary analysis of the acid mine water from a low 
grade uranium ore leaching factory indicates that the 
concentration of Mo (VI) ions is less than 100 mg/L; 



therefore in this paper, the effect of Mo (VI) ion concen- 
tration on the column performance was investigated in 
the range of 30 to 95 mg/L. The column was a simple 
glass tube with inner diameter of 1.5 cm and length of 
10 cm. Two plastic sieves both with pore size of 0.5 mm 
were installed at the top and bottom of this column. The 
experiments were conducted by pumping a metal solu- 
tion in up-flow mode through the packed-bed column 
with a peristaltic pump (Watson Marlow Pumps, Model 
205 U). A weight of 1 g of the dry biosorbent (equal to 
4 cm bed height) was packed within the column. 

Measurements and methods 

The functional groups of the biosorbent were determined 
by a Fourier Transform Inferred Spectrometer (Vector22- 
Bruker Company, Germany) in the range of 400-4000/cm 
wave number. The morphological analyses of the biosor- 
bent were characterized using a scanning electron micro- 
scope (SEM, JEOL JSM-6380). The specific surface area 
and pore volume were estimated using Brunauer- 
Emmett-Teller (BET) method. The analysis was carried 
out through an adsorption of N2 gas (gas effluent tem- 
perature of 300°C and bath temperature of 77.3°C). Micro 
pore volume was calculated from the volume of nitrogen 
adsorbed at relative pressure (P/Po) of 0.9893. Liquid 
samples were analyzed for residual metal ion concentra- 
tion by an inductively coupled plasma spectroscopy (ICP, 
Varian, Model Liberty 150 AX Turbo). Analytical wave- 
length was set at 202.031 nm for Mo (VI) ions. 

Analysis and mathematical modeling of breakthrough 
curves 

The design and optimization of packed-bed columns 
are difficult to carry out a priori without a quantita- 
tive modeling approach. From the perspective of 
process modeling, the dynamic behavior of a packed 
bed column is described in terms of the effluent 
concentration-time profile, i.e. the breakthrough 
curve. Dividing the metal mass (mad the quantity of 
metal retained in the column (mg)) by the sorbent 
mass (M) leads to the maximum uptake capacity of 
the biosorbent [25,26]. The breakthrough time, is 
said to be occurred when the effluent metal concen- 
tration reaches 5% of the influent value (hr) and bed 
exhaustion time, t^, is assumed to be occurred when 
the effluent metal concentration is equal to 95% of 
the influent value. Also, the effluent volume (Vg) that 
is the volume of metal solution passed into the col- 
umn (L) can be calculated as follows [25,27]: 

Ve = 0.06 -xQxte (1) 

where t^ is time (h) and Q is the volumetric flow rate 
(mL/min). Total amount metal sent to column m^tai 
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Figure 1 FTIR spectrum of (a) unloaded and (b) Mo (Vl)-loaded Ca-pretreated C. indica biomass. 



(mg) can be calculated as follows [25]: 



mtotal = 0.06 X CoxQxte 



(2) 



Total metal removal percent with respect to flow 
volume can be calculated as follows [25]: 



Total metal removal i 



yngd 

motal 



X 100 



(3) 



Successful design of a continuous packed bed column 
generally requires prediction of concentration-time pro- 
files for the effluents (i.e. breakthrough curves). The 
maximum uptake capacity of a biosorbent is also 
needed. There are several mathematical models in the 
literature, which have been used to represent the dynam- 
ics of the packed bed column. Thomas, Yoon- Nelson, 
Yan, and Belter models are among them, which can be 
written respectively as below [28-31]: 



Table 1 Physical properties of unloaded and Mo 
(Vl)-loaded Ca-pretreated C. indica biomass 



Property 


Unloaded biomass 


IVIo(VI)-loaded biomass 


Moisture (%) 


0.890 


0.740 


Density (g/cm) 


1.112 


1.241 


Pore volume (cm"^ g"^ 


5.259E-3 


2.700E-2 


Pore size (nm) 


2.871 


1.184 


Porosity (%) 


0.585 


3.350 


Specific area (mVg) 


36.710 


31.660 



C _ 1 

Co 1 + exp{KYNT - KyNt) 

C _ . 1 

Co" 



1 + 



C I 

Co~2 



1 + erf 



V2(T to. 



(4) 

(5) 
(6) 

(7) 



where Kxh is the Thomas rate constant (1/min. mg), qo is 
the maximum uptake capacity (mg/g). The kinetic coef- 
ficient kxh and the adsorption capacity of the bed qoTh 
can be determined from a plot of In ((Cq/C) -1) against 
^ at a given flow rate. 

M is the amount of adsorbent in the column (g), Cq is 
the influent metal ion concentration (mg L'^), C is the 
effluent metal ion concentration (mg/L), Q is the flow 
rate, Kyn is the Yoon-Nelson rate constant (1/min), t is 
the time required for 50% adsorbate breakthrough (min). 
The parameters kyN and x may be determined from a 
plot of In (C/(Co - C)) versus sampling time (t). t is the 
flow time (min) and "a" is the constant of the Yan 
model. 

Also, erf (x) is the error function of x, tohe temporal 
parameter (min) which indicates time needed for the 
outlet metal concentration to be the half of the one inlet 
metal concentration and a is the standard deviation of 
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Figure 2 SEM image of Ca-pertreated C. Indica biomass at different magnifications. 



the linear part of the breakthrough curve in the Belter 
model 

Results 

Adsorption mechanism of Mo (VI) on the C. indica 
brown algae 

Fourier Transform Infrared (FTIR) spectroscopy has 
proven and accepted to be a powerful tool for the study 
of biological samples [32]. Therefore, the algal samples 
were studied using FTIR technique. FTIR spectra of Mo 
(VI) unloaded and loaded forms of the biosorbent in the 
range of 400-4000/cm wave number were taken and 
presented in Figure 1. As shown in Figure 1, the 
unloaded biomass displays a number of absorption 
peaks, indicating the complex nature of the biomass 
examined. 



Also, some physical characteristics such as. Surface 
area, pore volume, and pore diameter of the Mo-loaded 
Ca-pretreated C. indica and unloaded Ca-pretreated 
C. indica biomass were determined on the basis of the 
Brunauer, Emmet and Teller (BET) method and the 
results were shown in Table 1. By the scanning electron 
microscopy (SEM), the biomass surface SEM image at 
different magnifications was obtained and is shown in 
Figure 2. 

The effect of flow rate 

The effect of flow rate on Mo (VI) biosorption by 
C. indica was studied by varying the flow rate from 
0.7 to 2.1 mL/min, while the bed height and initial 
Mo (VI) concentration were held constant at 4 cm 
and 30 mg/L, respectively. Figure 3 shows the 
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Figure 3 Breakthrough curves for Mo (VI) biosorption onto C. indica at different flow rates. 
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Table 2 Column data and parameters obtained at 30 mg/L inlet metal ion concentration and different flow rates for 
Mo (VI) biosorption onto Ca-pretreated C. indica biomass 

Q (mL/min) te (h) Vgff at te (I) tb (h) Uptake capacity (mg/g) Texp (h) Ad (%) 

0.7 17.14 0.72 0.06 10.32 8.16 47 

1.4 17.01 1.42 0.05 18.32 5.51 61 

2.1 9.05 1.14 0.02 11.7 2.34 35 



breakthrough profile of Mo (VI) biosorption at differ- 
ent flow rates and the results of the breakthrough 
curve analysis are given in Table 2. Results showed 
that the breakthrough and exhaustion time decrease 
with the increase of the flow rate from 0.06 to 
0.02 h and from 17.14 to 9.05 h respectively. Also, in 
this range, the metal uptake capacity was strongly 
influenced by increasing flow rate from 0.7 to 
1.4 mL/min, so that, its value for molybdenum 
increased from 10.32 to 18.32 mg/g. As can be seen 
from Table 2 when the flow rate was increased from 
1.4 to 2.1 mL/min' the molybdenum removal capacity 
was decreased from 18.32 to 11.7 mg/g. 



The effect of influent concentration 

The change in the influent metal ion concentration has 
a significant effect on breakthrough curve. The break- 
through curves obtained by changing inlet Mo (VI) con- 
centration from 30 to 95 mg/L at 4 cm bed height and 
1.4 mL/min flow rate are shown in Figure 4. The in- 
crease in the initial concentration of Mo (VI) solution 
increases the adsorption capacity from 18.32 to 
30.19 mg/g and decreases the percentage of Mo (VI) re- 
moval from 61 to 38% (Table 3). Also, the breakthrough 
time and exhaustion time decreased with an increase in 
initial concentration as shown in Table 3, the treated 



volume was the greatest (1.42 L) at the lowest inlet 
concentration. 

Model of column data 

The design and optimization of a fixed-bed sorption col- 
umn involves the employment of mathematical models, 
which must be used to describe and predict the 
experimental breakthrough curves, in order to make the 
scale up of the process possible. Because of this, the 
experimental breakthrough curves have been fitted to 
each one of the models mentioned previously with the 
aim to describe the behavior of the column for biosorp- 
tion of Mo (VI) onto C. indica and to determine the 
corresponding kinetics parameters. The experimental 
breakthrough curves have been fitted to the models 
through the non-linear regression, using the MATLAB 
software. 

Table 3 lists the model parameters of Thomas, Yoon 
and Nelson, Yan and Belter models with values. 
Figures 5, 6, 7 and 8 show the experimental points and 
predicted breakthrough curves in different conditions 
according to the Thomas, Yoon and Nelson, Yan and 
Belter models, respectively. 

As can be seen, these models (except the Yan model) 
were found to be most suitable to represent the kinetics 
of biosorption of Mo (VI) in a packed bed of C. indica. 
But the Belter model was better. 
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Figure 4 Breakthrough curves for Mo (VI) biosorption onto C Indica at different influent concentrations. 
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Table 3 Column data and parameters obtained at 1.4 mL/min and different Inlet metal Ion concentrations for Mo (VI) 
blosorptlon onto Ca-pretreated C. indica blomass 



Co (mg/L) 


te (h) 


Veff at te (1) 


tb (h) 


Uptake capacity (mg/g) 


Texp (h) 


Ad (%) 


30 


17.010 


1.420 


0.05 


18.320 


5.51 


61 


60 


10.570 


0.887 


0.04 


20.650 


4.42 


42 


95 


10.500 


0.882 


0.03 


30.190 


3.39 


38 



Discussion 

Adsorption mechanism of Mo (VI) on the C. indica brown 
algae 

Molybdenum exhibits multiple valences of +3, +4, +5 
and +6. The stable valence of +6 is important, not only 
for extraction but also for chemical adsorption mechan- 
ism. Predominating ionic species of Mo (VI), such as 
MoOi^' MoOl', HM0O4 and neutral species H2M0O4 
are reported to be present in a molybdenum aqueous 
solution [33]. Such species play a vital role to under- 
stand the biosorption mechanism of molybdenum onto 
C. indica biomass. Brown algae contain 20-40% of 
alginic acid on a dry weight basis. Alginic acid, a linear 



polysaccharide, is a main cell wall constituent of brown 
algae. Numerous chemical groups have been proposed 
to be responsible of biosorption metal binding by algae 
(hydroxyl, carboxyl, sulfhydryl, sulfonate, etc.); their im- 
portance for metal uptake depends on factors such as 
the quantity of sites, its accessibility, chemical state or 
affinity between site and metal [20]. The confirmation of 
the presence of these groups has been achieved by FTIR 
analysis. When the pH of solution decreases the hy- 
droxyl groups on the cell walls of biomass is protonated 
and the surface of biomass is positively charged that 
causes an increase in electrostatic adsorption of anionic 
species of molybdenum. A similar phenomenon has also 
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Figure 5 Comparison of the experimental and model fit breakthrough curves for Mo (VI) biosorption by C indica at different flow 
rates and influent concentrations according to the Thomas model. 



Kafshgari et al. Iranian Journal of Environmental Health Science & Engineering 2013, 10:14 
http://www.ijehse.conn/content/1 0/1 /1 4 



Page 7 of 1 1 




The Yoon-Nelson model C0=30 ppm 

The Yoon-Nelson model C0=60ppm 

The Yoon-Nelson model C0=95 ppm 

♦ Exp (C0=30 ppm) 
A Exp (C0=60 ppm) 
X Exp (C0=95 pp m) 



10 



20 30 
time (hr) 



40 



50 



60 




The Yoon-Nelson model Q=0.7 ml/min 

The Yoon-Nelson model Q=1.4 ml/min 

The Yoon-Nelson model Q=2.1 ml/min 

♦ Exp (Q=0. 7 ml/min) 
A Exp (Q=l. 4 ml/min) 
X Exp (Q=2.1 ml/min) 



10 



20 30 
time (111) 



40 



50 



60 



Figure 6 Comparison of the experimental and model fit breakthrough curves for Mo (VI) biosorption by C indica at different flow 
rates and influent concentrations according to the Yoon and Nelson model. 



been reported for Mo (VI) adsorption onto persimmon 
residua [34]. 

In order to further explain the adsorption mechanism 
of Mo (VI), FTIR spectroscopies of biomass before and 
after molybdenum adsorption were investigated. The 
FTIR spectra of unloaded biomass showed broad bands 
at 3476/cm that represent bounded hydroxyl (-OH) and 
amine (-NH2) groups. The bands observed at 2899/cm^ 
could are assigned to the aliphatic C-H group. The 
peak around at 2369/cm^represent stretching vibrations 
of -NH2, -NH^ and -NH groups of the biomass. The ad- 
sorption bands at 1656/cm and 1421 /cm were attributed 
to stretching vibration of carboxyl group (-C = O). The 
peak around at 1082/cm were due to the C-O stretching 
of alcohols and carboxylic acids. The bands at 893 cm"^ 
can be attributed to the aromatic C-H bending vibra- 
tions. Also, the strong band around 1100-1000/cm are 
due to the C-OH bond, which is a characteristic peak 
for polysaccharides [35]. After Mo (VI) biosorption, the 
stretching vibration bands of hydroxyl and amine groups 
were shifted to 3446/cm for Mo (VI) -loaded biomass. 
The vibration bands of aliphatic C-H group were shifted 



to 2972/cm. The stretching vibration bands of carboxyl 
groups were shifted to 1646 and 1377/cm. The bands 
assigned to C-O stretching and aromatic C-H bending 
vibrations were also shifted to 1007 and 879/cm, res- 
pectively. These changes observed in the spectrum indi- 
cated the possible involvement in Mo (VI) biosorption 
process of those ftinctional groups on the surface of the 
biomass. In addition, the disappearance of the band at 
2369/cm indicated that the chemical interactions between 
the protonated amine group of the biomass and molyb- 
denum ions was mainly involved in the biosorption. 

As can be seen from Table 2 the Mo (VI) biosorption 
process onto C. indica biomass led to change of density 
from 1.112 to 1.241 g cm"^. By comparison of biomass 
density with water density, we can observe that the 
biomass has a high degree of settle ability. The impli- 
cation is that the biomass would be suitable for the 
continuous flow systems. Also, specific area decreased 
from 36.71 to 31.66 m^/g and porosity increased from 
0.585 to 3.35%. These observations indicated that the 
chemical adsorption of Mo (VI) ions changes the cell 
wall structure of biomass. 
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Figure 7 Comparison of the experimental and model fit breakthrough curves for Mo (VI) biosorption by C indica at different flow 
rates and influent concentrations according to the Yan model. 



The effect of flow rate 

As can be seen from Table 2, the breakthrough curves 
became steeper and the breakthrough and exhaustion 
time decreased with the increasing flow rate. This 
behavior may be due to insufficient residence time of the 
solute in the column, which causes the Mo (VI) solution 
to leave the column before equilibrium occurs. On the 
other hand, the results showed that in the range of low 
flow rates (high residence times), the overall rate of Mo 
(VI) ions removal in the packed column was controlled 
by external mass transfer limitations. At low flow rates 
(from 0.7 to 1.4 mL/min), the metal uptake was strongly 
influenced by increasing the flow rate, so that, its value 
increased from 10.32 to 18.32 mg/g. In addition, the 
results showed that in the range of higher flow rates 
(lower residence times), the overall rate of metal sorp- 
tion by C. indica biomass was controlled by diffusion 
limitations of the solute into the pores of sorbent. When 
the flow rate increased from 1.4 to 2.1 mL/min, the 
liquid residence time in the column decreased, resulting 



in a lesser biosorption of metal ions, and hence, the 
metal uptake decreased from 18.32 to 11.7 mg/g. Similar 
observation was reported by Vijayaraghavan et al [36]. 

The sensitivity of the biosorption to the liquid flow 
rate can be explained by the fact that the liquid resi- 
dence time in the column is critical for the biosorption 
process. When the process is subjected to external mass 
transfer control, a higher flow rate decreases the liquid 
film resistance and when the process is subjected to 
intraparticle mass transfer control, a slower flow rate 
favors the sorption. 

The effect of influent concentration 

One of the parameters that strongly affect the metal 
removal capacity as well as the general position of the 
breakthrough curve is the initial metal ion concentration. 
According to the obtained results the uptake capacity 
increases as the inlet Mo (VI) concentration due to the in- 
creasing of driving force for mass transfer. This result is in 
agreement with that reported by Vijayaraghavan et al 
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Figure 8 Comparison of the experimental and model fit breakthrough curves for Mo (VI) biosorption by C indica at different flow rates and 
influent concentrations according to the Belter model. 



[36,37], Han et al [38] and Malkoc et al [39] studies. On 
the other hand, the percentage of Mo (VI) removal 
decreases as the influent concentration of Mo (VI) in- 
creases because, when the column gets close to saturation, 
the quantity of metal biosorbed is very low in comparison 
to the quantity of metal that passes through the column. 
A similar trend was also observed for Cr (III) biosorption 
by olive stone [40]. 

As influent concentration increased, the breakthrough 
curve became steeper and shifted towards the origin for 
Mo (VI) ions as the binding sites became more quickly 
saturated in the system at higher initial concentrations. 

Also, as shown in Table 3, the treated volume was the 
greatest at the lowest inlet concentration since the lower 
concentration gradient caused a slower transport due to 
a decreased mass transfer coefficient. 

Model of column data 

Thomas model has been applied to the breakthrough 
curve prediction. The calculated qo value is maximized 
in flow rate of 1.4 mL/min and 30 mg/g inlet metal 



Table 4 Thomas, Yoon-Nelson, Belter and Yan model 
parameters at different flow rates and inlet metal ion 
concentrations for Mo (VI) biosorption onto Ca-pretreated 
C. indica biomass 



Experimental 


Q (mL/min) 


0.7 


1.4 


2.1 


1.4 


1.4 


conditions 


Co (mg/L) 


30 


30 


30 


60 


95 


Thomas model 


qo 


9.910 


16.125 


9.571 


21.082 


23.678 






0.008 


0.020 


0.025 


0.007 


0.004 




r2 


0.979 


0.984 


0.973 


0.941 


0.975 


Yoon-Nelson model 


T(h) 


7.865 


6.399 


2.532 


4.183 


3.132 




Kyn 


0.266 


0.611 


0.756 


0.402 


0.391 




r2 


0.979 


0.984 


0.973 


0.941 


0.975 


Belter model 


0 


0.794 


1.188 


0.852 


0.987 


1.341 




to (h) 


9.895 


13.25 


3.020 


4.235 


2.349 




r2 


0.982 


0.990 


0.973 


0.943 


0.976 


Yan model 


Qo 


4.545 


15.4 


10.49 


13.7 


14.27 




a 


1.418 


1.903 


1.409 


1.087 


0.996 




r2 


0.956 


0.946 


0.946 


0.833 


0.881 
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concentration, which agrees with the one obtained by 
the experimental data (Table 4, Figure 5). 

According to the obtained results in Table 4, the sensi- 
tivity of the Thomas model constants with the changes 
of the liquid flow rates and the existence of a maximum 
in the qo values as the flow rate increases show that the 
controUed-rate step shifts from external to internal mass 
transfer limitations and the model is able to predict it. 
Also, the kinetics constant, Kxh? tends to increase as the 
flow rate increases and decrease as the influent metal 
concentration increases. 

For the Yoon and Nelson model, the calculated t 
values are quite close to those found experimentally, 
which indicates that the parameters of the model are 
similar to those obtained in the experiments (Tables 2, 3, 
4 and Figure 6). On the other hand, the time necessary 
to reach 50% of the retention, t, significantly decreases 
when the flow rate and influent metal concentration 
increase, because the saturation of the column occurs 
more rapidly. 

The equation proposed by Yan was not able to fit the 
experiment breakthrough curves in a non-linear analysis 
(Figure 7). The parameters obtained by non-linear fit of 
the data, qo, and a along with correlation coefficients is 
listed in Table 4. As can be seen from Table 4, the Yan 
model does not adequately reproduce the experimental 
data, obtaining values less than 0.956. Also, the qo 
values predicted by this model are lower than those 
found from the experimental data, especially in higher 
inlet molybdenum concentrations, which once again 
indicates that it is not possible to reproduce the results 
using this model. 

The Belter model parameters to.s and o for the 
biosorption of Mo (VI) on Ca-pretreated C. indica 
are found by fitting Eq. (7) to the experimental data 
(Table 4). Breakthrough curves calculated using these fit 
values for to.s and a are shown in Figure 8, in compari- 
son with the experimental data. The model showed a 
good fit to most of the obtained column data. However, 
the model parameters of the Belter show no apparent 
dependence on the flow rate and inlet metal ion concen- 
tration variables. It must be taken into account that 
these empirical model parameters lack any physical sig- 
nificance and are only valid under the same conditions 
applied to obtain the breakthrough curves, because they 
can change with variables such as the flow rate, inlet 
metal ion concentration, column dimensions, or packed 
algae. Nevertheless, this shortcoming is also present 
when the simplified mechanistic models are used. To 
summarize, it can be said that the Belter model is able 
to describe the complete breakthrough curve, suggesting 
that the features of the dynamics of sorption column 
have been taken into account by the mathematical form 
of the model equation. 



Conclusions 

This study indicated that the Ca-pretreated C. indica 
biomass could be used as an efficient biosorption bed in 
the fixed-bed columns for the removal of MO (VI) ion 
bearing aqueous solutions because, C. indica is an inex- 
pensive, easily available biomaterial with high adsorption 
capacity for MO (VI) ion. The effect of concentration 
and flow rate on biosorption process was investigated and 
the experimental breakthrough curves were obtained. 
Also, it was observed that there was an optimum flow rate 
of 1.4 ml min'^ for metal ion in the column. 

Thomas, Yoon-Nelson, Belter and Yan models were 
applied to the experimental data obtained from the bio- 
sorption of MO (VI) ion onto C. indica. Among these 
models, the Belter model appeared to describe the 
experimental results better. 

Thus, Cystoseira indica can be used in the fixed bed 
column as a potential biosorbent for treatment of Mo 
(VI) polluted aqueous solutions. 
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